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Abstract 

The completion of the third-order QCD corrections to the inclusive top-pair pro¬ 
duction cross section near threshold demonstrates that the strong dynamics is 
under control at the few percent level. In this paper we consider the effects of the 
Higgs boson on the cross section and, for the hrst time, combine the third-order 
QCD result with the third-order P-wave, the leading QED and the leading non¬ 
resonant contributions. We study the size of the different effects and investigate 
the sensitivity of the cross section to variations of the top-quark Yukawa coupling 
due to possible new physics effects. 


1 Introduction 


Top anti-top quark production near threshold in collisions provides a unique op¬ 

portunity to measure the top-quark mass precisely, due to the well-dehned center-of-mass 
energy and the enhancement of the cross section due to the strong-interaction Coulomb 
force. Whether the required theoretical precision on the cross section can be achieved has 
been an open question, since the second order (non-relativistic) QCD calculations revealed 
unexpectedly large corrections and uncertainties inEi. After many years of work, the 
third-order QCD calculation has been recently hnished |^, resulting in a largely reduced 
theoretical uncertainty. With QCD effects under control, the emphasis shifts to other ef¬ 
fects which must be addressed for a realistic cross section prediction. The most important 
are Higgs effects associated with the top-quark Yukawa coupling, general electromagnetic 
and electroweak corrections, non-resonant production of the hnal state W^W~bb in the 
center-of-mass region near twice the top-quark mass, 2mt, and photon initial-state radia¬ 
tion. 

In this paper, we mainly focus on Higgs-boson effects and the sensitivity to the top- 
quark Yukawa coupling. The Yukawa potential generated by Higgs exchange |1] and one- 
loop corrections to tt production [5] have been considered long ago, but these early cal¬ 
culations do not reach the precision that corresponds to the third-order QCD calculation 
in the non-relativistic power-counting scheme. Third-order Higgs corrections to the pro¬ 
duction vertex and the energy and wave-function at the origin of a hypothetical S-wave 
toponium resonance have been computed in [S], but the ti cross section has not yet been 
considered. We supply this missing piece here. We also add for the hrst time the P-wave 
contributions [7] and the leading non-resonant contributions [8l|9] to the third-order S-wave 
QCD calculation. We then allow the top-quark Yukawa coupling yt to deviate from the 
Standard Model relationEl rrit = ytv/\/2 and investigate the sensitivity to such deviations 
given the current theoretical uncertainties. 


2 Higgs effects at NNNLO 

The contribution of the Higgs boson to the top pair production cross section e''“e“ —)■ ti 
introduces two new parameters, the Higgs mass mn, and top-quark Yukawa coupling yt. 
To set up the calculation we have to £x their relation to mt and the strong and electroweak 
couplings, as and ctew; to establish the power counting. Recall that it is customary to 
count as ^ V and aEW ~ ctgj where v = [(\/s —is the small top-quark velocity. 
A contribution of order (or, equivalently, v^) according to this counting is called “N*^LO” 
or “/cth order”. We opt for counting y^ ~ oew ~ and uih ~ Other options would 
be to count the top-Yukawa coupling like the strong coupling, y^ ^ as, or the Higgs mass 
ruH ~ fntv, or both. Clearly, with mt ~ 173 GeV, mn ~ 125 GeV and v ~ 1/10, the 
counting mn ~ is more appropriate. In the terminology of non-relativistic effective 

* The symbol v is used for the Higgs vacuum expectation value and the top-quark velocity, see below. 
The meaning should be clear from the context. 
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Figure 1: One-loop Higgs correction to the colour Coulomb potential. 


theory and the threshold expansion, the Higgs mass is of order of the hard scale, and 
not the potential scale, which has significant impact on the structure of the contributions. 
On the other hand the counting of the coupling simply determines at which orders in the 
expansion the Higgs contributions appear and we will justify our choice below. 

The effective held theory setup is described in detail in |T0]. We recall that the domi¬ 
nant S-wave production cross section is proportional to the imaginary part of the spectral 
function of the vector current 

= ^clG{,E) ^ , ( 2 . 1 ) 


where c„ is the hard matching coefficient of the vector current, E = — 2mt, and GiE) 

is the Green function in potential-nonrelativistic QCD (PNRQCD), i.e. the propagator 
of a non-relativistic top anti-top pair. The Higgs contributions to c„ are discussed in 
Section I2.1I To compute the corrections to the Green function the Higgs contributions 
to the PNRQGD Lagrangian have to be determined. Gounting mu ~ mt implies that 
the Yukawa-potential exp(—mj:^r)/r generated by Higgs exchange between the top quarks 
is replaced by the local interaction {y) j as is apparent from the Higgs propagator 
l/(q2 -|- in momentum space, where ~ can be neglected (expanded) relative 
to vEii- On the other hand, with mu ~ "nitVi both terms would have to be kept. The 
contribution to the momentum-space potential is therefore simply 


bnV 



( 2 . 2 ) 


We note that this is suppressed by with respect to the leading QGD Goulomb potential 
cts/q^, where one power of v arises from the counting of the Yukawa coupling, and two 
powers from the relative factor /m?^. The Higgs-induced potential is thus a NNNLO 
effect. The corresponding correction to the Green function G{E) is computed in Section lT^ 
Furthermore we have to consider corrections to the colour Goulomb potential as shown 
in Fig. [H With mu ~ mt counting, only the hard loop momentum region can yield a 
contribution. Since the external momenta are potential they have to be expanded, and we 
are left with an 0{yf) zero-momentum transfer correction to the top-quark-gluon 

coupling of the NRQGD Lagrangian. However, since the top field is renormalized in the 
on-shell scheme this contribution cancels. 
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2.1 Short-distance effects 


With yt ~ Os the leading, one-loop Higgs contribution to the hard matching coefficient 
of the vector current is of second order. It has been computed in [5l|Tll[T2]. Due to an 
additional diagram involving a. ZZH vertex, the matching coefficients for the jit and the 
vector part of the Ztt vertex differ. We neglect the contribution from this diagram and 
use the 'jtt coefficient, since the difference amounts to less than one percent of the already 
small NNLO Higgs contribution to the production vertex [5]. At NNNLO there are mixed 
Higgs and QCD corrections to the vector current. They have been computed as expansions 
for niH ~ fTit or mn 3> mt in [6], which is consistent with the adopted hard mn ~ 
scaling. We use the result expanded in (1 — which was denoted 16 in [6] and was 

shown to converge quickly for Higgs boson masses around the physical value of about 125 
GeV. Based on the results of [B] we estimate the truncation error due to this expansion 
to be well below one percent of the NNNLO matching coefficient and neglect it in the 
following. The NNNLO correction to the hard matching contains an IR divergence, which 
can be absorbed into a renormalization constant for the vector current 

- 1 + [pure QCD] + —(2.3) 


The renormalized hard matching coefficient for the vector current can be parametrized by 


= 1 -|- [pure QCD] -|- 
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(2.4) 


where can be obtained from [6] . For convenience we reproduce the relevant expressions 
in the Appendix. The pure QCD correction is also known to NNNLO [13]. Starting at 
NNLO there is also a self-energy correction to the Z boson mediated cross section that 
contains the Higgs boson. Since it does not involve the top Yukawa coupling we do not 
consider it here. It will be added together with other NNLO electroweak and non-resonant 
effects in future work. 

To justify the power counting yt ~ as we compare the size of the Higgs effects discussed 
above to their QCD counterparts. For the hard matching coefficient we obtain, with 
as{y = 80 GeV) = 0.1209, 


= 1 - 0.103U^ - 0.022|„2 + 0.031|j, 2 - O.OTO^s - 0.019|y2„^ + ..., (2.5) 

where the contributions from different orders of the couplings are shown explicitly. The 
power counting is clearly valid here. For the potentials it is natural to compare the Higgs 
potential 02.21) to the spin-projected QCD NNLO Darwin potential 5Vu, which is also 
local. Adopting as from above, we hnd 6hV = —0.98/m^, which is only slightly smaller 
than 5Vd = SnasCp/i^mf) = 1.35/m^. However, since the Darwin potential yields only a 
small correction compared to other NNLO effects we conclude that the overall counting is 
appropriate. 
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2.2 Potential contributions 


The potential correction to the Green fnnction can be obtained by qnantnm mechanical 
(PNRQCD) pertnrbation theory due to the instantaneous, hence particle number conserv¬ 
ing nature of potentials. Since the Higgs potential fl2.2p is a NNNLO effect, only the single 
insertion of 6hV is required to compute the NNNLO correction to the Green function 

6hG{E) = {0\Go{E)z6HVtGomO) = -SnVGoiEf. (2.6) 


The remarkably simple form arises because of the locality of the potential. The Green 
function Go{E) describes the propagation of a top quark pair, produced and destroyed at 
zero spatial separation, under the influence of the leading-order QGD Goulomb potential. 
The insertion of a local interaction thus factorizes into the product of a Green function 
to the left and the right of the insertion. Using the well-known result for the LO Green 
function 


Go{E) 


mfasGp 

Att 


1^11 


<j(l-A) + 0 (e) , 


(2.7) 


in d = 4 — 2e space-time dimensions, expressed through A = agGp/i^^—E/nrit), Lx = 
log^Xfi/^mtasGp)), and = 7 e + we observe that the imaginary part of (12.61) is 
UV divergent, 

Im |feG(i5)||,i, = Im |G„(B)], (2.8) 

where Go{E) denotes the exact d-dimensional LO Green function. Note that this expression 
also contains a finite term from the 1/e pole multiplying the unknown 0(e) part of the LO 
Green function. However, (12.8p exactly cancels in the combination Im[(c^Z“^)^G(U)]. This 
is analogous to the cancellation between the divergent part of the QGD Darwin potential 
single insertion and part of the divergence of the two-loop QGD contribution to c^. In the 
following, we therefore only have to consider the finite part 


^nGf^niE) 


yt 

2m\j 


m^asGp 

dvr 


n + i-T-v)(i-A) 


(2.9) 


Due to the non-perturbative treatment of the LO Goulomb potential in PNRQGD, 
the exact Green function contains single poles below threshold, which correspond to ^Si 
toponium bound states: 


G{E) 


E—^En 


En — E — ie 


( 2 . 10 ) 


The energy levels En and squared wave functions at the origin, |'0„(O)p, also receive cor¬ 
rections from the insertion of the Higgs potential. We use the parametrization 


E„ = Ei"' L + [pure QCD] + ' 

l^n(0)P = + [pure QGD] + , 


( 2 . 11 ) 
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where 





( 2 . 12 ) 


The corrections e//, fn can be obtained by expanding (I2.9p and (I2.10p around the bound- 
state energies and comparing coefficients. Alternatively both equations can be expanded 
around positive integer values of A. For fl2.9l) we obtain 



where = log{nfi/{mtasCp)) and the ellipsis denotes terms that are regular in the limit 
\ ^ n. We obtain 




(2.14) 


Gh = - 2 - 

mjj n 


Here Si (n) = ^ ^ denotes the harmonic number of order one. The result for en 


agrees with [6] and for fu we reproduce the value for n = 1 given in [6]. 

2.3 Combined 

Combining hard and potential effects due to the Higgs boson, more precisely, the top- 
Yukawa interaction with the Higgs boson, we can express the NNLO correction to the 
vector correlation function (12. ip as 



(2.15) 


The NNNLO correction is 


^ c®Ci)Go(£) + ^,SS^G(E) + 6„GUE)] . (2.16) 


It includes cross terms of the NNLO Higgs correction with the known NLO QCD corrections 
to the Green function, G{E) = Gq{E) -|- 5iG{E) and matching coefficient c„ = 

1 -t- Q;sCi/(47r) -|- .... Both terms are hnite, as it is understood here that Gq{E) from fl2.7p 
is minimally subtracted. The physical cross section is related to the imaginary part of 




3 Non-resonant and QED effects 


In the computation of QCD and Higgs corrections to the tt cross section the top decay 


width has been accounted for by the replacement E —)■ E + iTt, where both quantities are of 


the same order. Since the top quark is unstable, one should rather consider the production 
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cross section for the decay product W^W~bb of the top pair|i| This dual state can also be 
produced without an intermediate resonant top pair and only the sum of both processes 
constitutes a physical quantity. This is also apparent from an incomplete cancellation of 
UV divergences and scale dependence in the non-relativistic description of the resonant 
process starting at NNLO. The non-resonant correction is important for realistic cross 
section predictions, since it affects particularly the cross section below the position of the 
peak, where the sensitivity to the top-quark mass is largest [3]. 

The computation of the e’''e“ —)■ W^W~bb cross section in the top anti-top threshold 
region, such that it is consistent with the non-relativistic power counting, expansion and 
resummation of the resonant sub-process, can be performed in the framework of unstable- 
particle effective theory [TU[T5]. The master formula for the cross section is 

tA = f d^x(e-e+|T[zC>f)t(0)*C)W(a^)]|e”e+) 

k,l 

+ ‘^4?(e"e+|iC>ig^(0)|e"e+). (3.1) 

k 

The hrst line describes the resonant production of the tt pair through an operator Op^\ 
Since the initial state is colour-neutral, the matrix element further factorizes into a leptonic 
and a hadronic tensor as long as only QCD and no electroweak effects are considered, and 
the non-relativistic treatment in terms of the PNRQCD Green function G{E) is recovered. 
The second line describes non-resonant production of the W^W~bb hnal state, which 
is why the hadronic contribution can be absorbed fully into a hard Wilson coefficient. 
The leading non-resonant effects appear already at NLO in the non-relativistic power 
counting, and have been determined in |8], including the possibility of imposing invariant- 
mass cuts on the top decay products. The result for the total cross section was confirmed 
in [9] with an independent method. In the following analysis we combine the results 
from [S] with the third-order QCD calculation [5]. We exclude the small contribution 
from —)■ W^W~H followed by if —)■ bh, since it can be considered as a reducible 

“background” and eliminated by an invariant-mass cut on the bb jets as discussed in [8]. 
At NNLO only partial results for the non-resonant term in the second line of fl3.ip are 
available We do not consider them here and hope to include the complete NNLO 

non-resonant correction together with other NNLO electroweak effects in future work. We 
note that the cancellation of divergences at NNLO has already been demonstrated [13] and 
that the contribution is likely to be numerically relevant below threshold. 

There is, however, one further electroweak effect already at NLO. The counting a ~ 
of the QED coupling implies that the QED Coulomb potential (5Vqed = —47raQ^ /repre¬ 
sents a NLO correction relative to the leading order QCD Coulomb potential. While only 
a single insertion of this potential would be required for NLO accuracy, we include con¬ 
tributions involving the NLO QED Coulomb potential up to NNNLO, i.e. we also include 
multiple insertions of 5 Vqed as well as mixed insertions together with other potentials and 
current matching coefficients. The required expressions can be inferred from the known 

Whe W boson can be treated as stable here, since its kinematics is not sensitive to the W width. 
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results for multiple insertions of the QCD potentials [HI[19], but one has to be careful 
when considering insertions which contain divergences in the imaginary part, since the 
QCD Coulomb potential contains 0{e) terms which are absent in the QED potential. The 
QED potential has already been included in previous calculations [8l|2ni[2T], but not yet 
in combination with the third-order QCD result [3]. Similarly, the numerical effect of the 
NLO non-resonant terms was studied in detail in [8], but was not implemented so far in 
the code that includes the higher-order QCD corrections. 

4 Size of Higgs and other non-QCD effects 

For the cross section predictions shown below we employ the values 

mf®(20 GeV) = 171.5 GeV, E* = 1.33 GeV (4.1) 

for the top-quark PS mass [22] and top-quark width. We note that the QED Coulomb 
potential is not part of the definition of the PS mass, since, in contrast to QCD, higher- 
order QED corrections do not give rise to an IR renormalon ambiguity in the pole mass, 
and are rapidly decreasing. The strong and electromagnetic couplings are 

as{Mz) = 0.1185 ± 0.006, a{Mz) = 1/128.944, (4.2) 

where the QCD coupling refers to the MS scheme and the running QED coupling, taken 
from [23], to the on-shell scheme. These parameters are taken to be consistent with [3]. 
We further use 


Mw = 80.385 GeV, Mz = 91.1876 GeV, mn = 125 GeV (4.3) 

for the electroweak gauge and Higgs boson masses, from which we derive the Weinberg 
angle, Higgs expectation value v and top-Yukawa coupling yt = \/2mt/v through tree- 
level relations. Here rrit is the top pole mass, computed from the PS mass with NNNLO 
accuracy. 

The QGD NNNLO result includes the P-wave contribution at the same order, computed 
in [7], which arises from production of the tt pair through the axial-vector coupling of a 
virtual Z-boson. This enhances the S'-wave cross section presented in [5] by about 1%. The 
theoretical uncertainty of the cross section calculation itself is estimated by varying the 
renormalization scale between 50 and 350 GeV. The “default scale” is set to 80 GeV. The 
“finite-width” factorization scale related to the separation of resonant and non-resonant 
terms in fl3.1l) is hxed to = 350 GeV. The dependence on this scale is cancelled exactly 
order-by-order in the sum of the two contributions. Since presently the non-resonant terms 
are included only to NLO, while the resonant terms are known to NNNLO, there is a small 
uncancelled dependence on at NNLO. 

In Fig. |2]we show the total e+e“ —)■ W^W~bb cross section in the range of e+e“ center- 
of-mass energy ^/s a few GeV below and above the top anti-top threshold including the 
Higgs, QED and non-resonant corrections discussed above (red/dark-grey hatched band) 
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Figure 2: The R ratio for QCD-only (light-grey hatched band) and including the Higgs, 
QED, and non-resonant contributions (red/dark-grey hatched band) as functions of the 
center-of-mass energy. The bands are due to the scale variation. The upper plot shows the 
absolute results and the lower plot the results normalized to the full one evaluated at the 
scale /i = 80 GeV. 
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Figure 3: The sum of the QCD and individual non-QCD contributions to the R ratio, 
normalized to the QCD-only result, as functions of the center-of-mass energy. 


and compare it to the QCD-only result at NNNLO (light-grey hatched band). The cross 
section is always normalized to the LO cross section for e’''e“ —)■ yielding the so- 

called R ratio R = a{e~^e~ —)■ W~^W~bb)/ao, with do = 47ra^/(3s). The QED and Higgs 
potentials lead to an attractive force, which enhances the cross section. The leading Higgs 
contribution to the short-distance coefficient (12.hh is also positive, resulting in an overall 
enhancement of about 10% near and 5% above the peak. Below the peak the negative non¬ 
resonant contribution becomes important and wins over the QED and Higgs enhancement. 
The lower plot in Fig. |2] contains the same results as the upper one, but now all values have 
been normalized to the full R ratio evaluated for /x = 80 GeV. Comparing the two bands, 
we again observe the enhancement of the full cross section around and above the peak. 
The rise in the QCD-only result at energies below the peak is due to the non-resonant 
contribution, which decreases the full cross section in this region. We observe a small 
increase in the scale uncertainty of the full result relative to QCD-only, which can now 
reach ±5% about one GeV below the peak, but is mostly of the ±3% size. The additional 
scale dependence arises mainly from the Higgs potential insertion. 

The size of the three non-QCD contributions considered here— the Higgs, QED and 
non-resonant contributions— are shown separately in Fig. [3l normalized to the NNNLO 
QCD-only result for the total ti cross section and plotted as function of the center- 
of-mass energy s/s (always fi = 80 GeV). The peak of the QGD-only cross section for the 
adopted parameters is at y/s = 343.95 GeV (see Fig.Hj). The Higgs and QED contributions 
result in positive shifts of the cross section of about 4-8% for the former and 2-8% for 
the latter, depending on the value of i/i. Both shift the peak to lower energies, thus 
resulting in larger enhancements below the peak of the QGD-only result. The NLO non¬ 
resonant correction is a nearly energy-independent, negative contribution [8] in absolute 
size and is therefore increasingly important below threshold, which explains the shape of 
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the corresponding line in Fig. [3l Its absolute size is smaller than the sum of the QED and 
Higgs contributions in the peak region and above, which leads to the overall positive shift in 
these regions observed in Fig.[2J Below threshold the resonant contribution falls off quickly 
and the relative correction from the non-resonant part becomes very large, up to 20%. The 
same behaviour is found for the dependence on the scale /x^. For variations in physically 
reasonable ranges from the potential to the hard scale we hnd a relative uncertainty of less 
than ±1% above threshold, but up to ±3% a few GeV below threshold. The dominant 
part of this uncertainty comes from the NNLO corrections and cancels exactly once the 
NNLO non-resonant corrections are included, after which the remaining /x^ dependence 
is at most 1%. We note the uncertainty from /x^, variation is not included in Fig. [2] and 
subsequent figures that show scale variations. 

For a determination of the top-quark mass from the threshold cross section, the sharp 
rise and the peak in the cross section are its most important features. On the other hand, 
determinations of the top-quark decay width and its Yukawa coupling require a precise 
knowledge of the overall normalization. In order to judge the influence of the Higgs, QED, 
and non-resonant contribution on these quantities. Fig. 0] shows the impact on the peak 
(upper plot) and maximal slope (lower plot) of the cross section, when these contributions 
are added successively to the QCD-only result. The theoretical uncertainty due to the 
variation of as{Mz) within its uncertainty given in fl4.2p is shown as the inner error bars. 
The outer error bars are the quadratic sum of the as and scale uncertainty. They provide 
an indication of the significance of the changes. 

The Higgs and QED contributions result in a negative shift of the peak position and 
an increase in the peak height. Correspondingly, the position of the maximal slope is also 
shifted to a lower energy and its value is increased when these contributions are added to 
the QCD result. The Higgs contribution shifts the peak position by —35 MeV and the 
QED contributions adds another —71 MeV. Since the peak position is related to twice 
the top-quark mass, this translates into a —53 MeV difference in the top-quark masses 
obtained from the full and QCD-only results for the cross section. Since the non-resonant 
contribution is an almost energy-independent negative shift, it has almost no influence on 
the position of the peak and only decreases its height. It also leaves the slope unchanged. 
Therefore, it is mostly important for the overall normalization of the cross section. 

5 Sensitivity to the top Yukawa coupling 

The mechanism of fermion mass generation in the Standard Model (SM) is intimately 
related to the question whether the Yukawa coupling of the Higgs boson to fermion / is 
proportional to the fermion’s mass, yj = \/2mf/v. In the SM effective theory including 
dimension-six operators |2T1|25] this relation can be violated, for example, by the operator 

A/: = -^(0V)(Q30^ix) + h.c., (5.1) 

where cnp is a new, independent coupling, A the scale of new physics and 0 = For 

simplicity we have neglected flavor indices and assume the new physics to only affect the 
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Figure 4: The upper plot shows the peak height versus the peak position of the R ratio 
for QCD-only (red), QCD+Higgs (green), QCD+Higgs+QED (blue) and including all 
contributions (black). The inner error bar denotes the uncertainty due to as{Mz), while 
the outer one denotes the quadratic sum of the scale and as uncertainties. The lower plot 
is the same but for the maximal slope. In this case there is no visible change due to the 
non-resonant contribution. 
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third generation. After spontaneous symmetry breaking, the operator generates corrections 
to the top mass term and Higgs coupling, 

2 

D - + h.c., (5.2) 

where h denotes the physical Higgs held. We observe that the coefficients of the mass and 
Yukawa term differ, and obtain the relation 

K = = 1 + Si 

where 

( SM . 1 SM I SCNPU^ (r A\ 

Below we will use Kt dehned in fl5.3p to parametrize corrections to the standard relation 
between yt and m*. To investigate the sensitivity of the top anti-top cross section to Kf, we 
do not use the SM relation between rrit and yt in the calculation of the Higgs contribution, 
and rescale yt by Kt in the Higgs potential (12.2p and the short-distance contributions. 
That is, we simply assume that some new physics effect makes the top mass and Yukawa 
coupling independent parameters. Evidently, the complete set of dimension-six operators 
may induce further anomalous couplings of the top quark, such as an anomalous top-gluon 
coupling, which can give additional short-distance and potential contributions to the cross 
section. A full treatment of these effects is beyond the scope of this work. 

The sensitivity of the R ratio to variations of the Yukawa coupling is shown in the upper 
plot of Fig. [5l The plot shows curves for different values of Kt normalized to the result at 
Kt = 1 and /i = 80 GeV. The main effect of an increase (decrease) in the Yukawa coupling 
is a strengthening (weakening) of the attractive potential between the top and anti-top 
quarks. This results in an increase of the cross section of 5-10% for Kt = 1.5 or a decrease 
of 3-5% for Kt = 0.5, with some dependence on the center-of-mass energy. In order to 
provide a first estimate of the possible precision of a Yukawa coupling measurement from 
top anti-top threshold production, the plot also shows the theoretical uncertainty due to 
the variation of the renormalization scale y. Naively one would expect to be only sensitive 
to values of the Yukawa coupling that lie outside this uncertainty band. From the figure 
we see that this requires rather large deviations from the SM value of roughly -|-20% or 
—50%. However, a more detailed analysis should also take into account the shape of the 
curve, which may lead to an improved sensitivity. 

Another important point is that a variation of the strong coupling leads to similar 
changes in the cross section as a variation of the Yukawa coupling. This can be seen 
by comparing the upper and lower plot in Fig. [5l where the lower one shows curves for 
as{Mz) = 0.1205 and as{Mz) = 0.1165. Just as for the Yukawa coupling, an increase 
(decrease) of the strong coupling leads to an increased (decreased) cross section, though 
the energy dependence of the shift is slightly different. To make this point clearer. Fig. [6] 
shows the change in height and position of the peak of the R ratio due to changes in 
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Figure 5: The sensitivity of the R ratio to the variation of the Yukawa coupling (upper 
plot) and the strong coupling (lower plot). The bands denote the uncertainty due to scale 
variation of the R ratio for yt = yf^ (k^ = 1) and as{Mz) = 0.1185. All values have been 
normalized to the R ratio evaluated at /i = 80 GeV. 
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Figure 6: Changes in peak height and position due to variation of the Yukawa coupling 
(red line) and the strong coupling (green line). The black error bars denote the and 
combined scale and Ug uncertainty for yt = yf^ (k^ = 1) and as{Mz) = 0.1185 (cf. upper 
plot in Fig. H]). 


the strong and Yukawa coupling. The similar slope of the resulting lines indicates the 
degeneracy in the variations of the two parameters. Thus, the precision of a Yukawa 
coupling measurement depends on the uncertainty of as{Mz)- Alternatively, one could 
perform a simultaneous fit of both couplings, though this would again lead to a loss in 
precision for the Yukawa coupling. 

Ref. [26j finds that for a Higgs boson with a mass of about 125 GeV the top quark 
Yukawa coupling can be obtained with a statistical uncertainty of only 4.2%. This result 
is based on the increase of the cross section from yt = 0 to yt = yf^, which is assumed to 
be 9% and energy independent. Neither the theoretical uncertainty of the cross section, 
nor the correlation with the strong coupling constant are considered. Our results show 
that once theoretical uncertainties are taken into account, it is unlikely that such a high 
precision can be achieved. 


6 Conclusion 

The completion of the NNNLO QCD correction to the top anti-top production cross sec¬ 
tion near threshold has increased the precision of the theoretical prediction to a level where 
non-QCD effects gain importance. In this paper we added NNNLO Higgs, and the leading 
(NLO) QED and non-resonant contributions to the third-order QCD result for the top 
anti-top production cross section. All three effects are larger than the current QCD uncer¬ 
tainty of about ±3% [3] and cause a distinct modification of the cross section below, near 
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and above the peak. We quantified the theoretical uncertainty in the presence of these 
effects, the dependence on the strong coupling, and the sensitivity to a modification of the 
top-Yukawa coupling. Further studies should be performed in the framework of realistic 
simulations accounting for beam and initial-state radiation effects. On the theoretical side, 
the inclusion of NNLO electroweak and non-resonant corrections would further sharpen 
the prediction, especially below the resonance peak. 
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A Higgs contribution to the hard matching coefficient 

For convenience we give the expressions used for the Higgs contribution to the hard match¬ 
ing coefficient of the vector current fl2.4p . The results are taken from [6], only the prefactors 
have been adjusted to match our convention. 
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where y = 1 — z. 
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